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Abstract 

Transport properties in gases are significantly affected by temperature. In previous works it has been 
shown that when the thermal agitation in a gas is high enough, such that relativistic effects become 
relevant, heat dissipation is driven not solely by a temperature gradient but also by other vector forces 
m-m In the case of relativistic charged fluids, a heat flux is driven by an electrostatic field even in the 
single species case |8]. The present work generalizes such result by considering also a magnetic field in an 
arbitrary inertial reference frame. The corresponding constitutive equation is explicitly obtained showing 
that both electric and magnetic forces contribute to thermal dissipation. This result may lead to relevant 
effects in plasma dynamics. 


1 Introduction 

Linear irreversible thermodynamics predicts linear relations between dissipative fluxes and thermodynamic 
forces. These constitutive, or closure, equations when generalized for high temperature fluids feature cross-like 
effects. In particular, it has been shown that a heat flux is driven by an external electrostatic field in a charged 
single component fluid. This corresponds to an electro-thermal effect in simple gases which is not present in 
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the non-relativistic, low temperature, case when the full Boltzmann equation is considered |2]. In the present 
paper such an effect is generalized to the electromagnetic case considering an additional magnetic field and 
working in an arbitrary reference frame. 

Boltzmann’s equation has been applied by other authors to address both the single component charged 
gas as well as the mixture in both relativistic and non-relativistic scenarios. In particular, Spitzer established 
a heat flux depending on the electrostatic field for a monocomponent fluid by neglecting the time derivative 
term, assuming a steady state situation [a [IQ]. On the other hand, non-relativistic two component plasmas 
have been addressed widely in the literature where electro-thermal effects correspond to cross effects that rely 
on the system being a mixture Kremer et. ah adressed the relativistic charged mixture obtaining 

Fourier and Ohm’s law as well as electro-thermal effects driven by the electric field 

The formalism here presented follows the standard steps in the Chapman-Enskog method whithin Marie’s 
relaxation approximation for a charged single relativistic fluid. Once the electromagnetic contribution to the 
non-equilibrium part of the distribution is established, the heat flux is calculated based on the expression in 
terms of the energy flux arising solely from the chaotic component of the velocities na. The result obtained 
consists on a general expression for the heat flux present in a simple fluid driven by an external electromagnetic 
field in an arbitrary reference frame. 

The rest of the paper is divided as follows. Section II is devoted to the description of the fundamental 
elements of relativistic kinetic theory where Boltzmann’s equation is written for a charged single component 
fluid and the corresponding transport equations are obtained. A simplified collision model is incorporated 
in Sect. Ill in order to establish the basic structure of the non-equilibrium part of the distribution function. 
The heat flux is obtained in Sect. IV, where the explicit constitutive equation is shown and the corresponding 
transport coefficient calculated. The explicit dependence of the heat flux on the electromagnetic force is also 
shown. Conclusions and final remarks are included in Sect. V. 


2 Transport equations for a relativistic charged simple fluid 

Relativistic molecular dynamics become relevant in gases when thermal agitation is such that the particle’s 
chaotic velocities become close to the speed of light. The relativistic nature of the gas is measured through 
the relativistic parameter z = kT/me?. When z >1, the thermal energy is comparable to the rest mass of the 
individual molecules and relativistic corrections to particle dynamics permeate, through statistical averages, 
to measurable corrections in state variables and their evolution equations. 
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From the kinetic theory point of view, the formal treatment of high temperature gases is based on the 
relativistic Boltzmann equation. For the case of a single charged relativistic fluid, such equation can be 
written as PE] 




(1) 


dx°^ m dv 

Here / is the relativistic distribution function which, as in the non relativistic case, yields the occupation 

number of phase-space cells with volume d^xd^v. As long as the fluid is dilute and no intense gravitational 

fields are present, a flat Minkowski metric (with a I— signature) is appropriate where = [r, ct] and 

/ £ \-V2 

= 7 ^ [w^,c], with 'jw = (1 ~ ^ 7 ^) being the Lorentz factor. Here and in the rest of the work the 
Einstein summation convention is used, with greek indices running from 1 to 4 and latin ones up to 3. Also, 
in Eq. ([T]) the external force is expressed in terms of the electromagnetic field tensor Notice that in the 
electromagnetic case, the inclusion of a velocity dependent force field does not alter the streaming side of the 
equation because of the antisymmetry properties of that is 


d 

dv^ 




,dv^ 

dv° 


= 0 


( 2 ) 


Conserved quantities and the expressions for dissipative fluxes and state variables are obtained following 
a similar procedure as the one carried out in Ref. na. Starting from Boltzmann’s equation (Eq. ([T])) one 
multiplies both sides by a collisional invariant "0 (n^) and integrate over velocity space, which yields 


d 

dx^ 


'ip {v^) v°‘fd*v 


(n^) 



J iff) d*v 


(3) 


where the fact that velocity and position are independent variables in phase space has been used for the first 
term on the left hand side. The invariant volume element is given by d*v = 'y~^d^v [T7]. For the second term 
on the left hand side of Eq. ([3]) one has 


dv°‘ 


dv 




dv‘^ 


]d*v- [ fF^’'v^,^[lpivf]d*v 


dv^ 


(4) 


where the first term vanishes using Gauss’ theorem and the fact that the distribition function vanishes 
(exponentially) on the boundaries of the v domain. The second term also vanishes because of the antisymmetry 
of the field tensor. For the right hand side, the symmetry properties of the collision operator, together with 













the fact that V’ (v^) is a collisional invariant, leads to the vanishing of snch term and thus one obtains 


dx 


^ ^ iP (n“) v^fd*v - — [ = 0 


m 




( 5 ) 


Considering ip (v^) = 1 in Eq. ([S]) one has 


dN^ 

dx^ 


= 0 


( 6 ) 


with = J v°‘fd*v being the particle four flux. The relation of this flux with the hydrodynamic four velocity 
W" follows as in the neutral case as 


= nW 


(7) 


and thus, the conservation of such quantity leads to the continuity equation namely. 


dn dhi^ 

— 

dx^^ dxf^ 


( 8 ) 


For the momentum-energy invariant ^p = mv^ one obtains the following balance equation 


dT^c 

dx^ 


= qrdA^F^'' 


(9) 


where the fluid’s energy-momentum tensor is defined as 


= / mv^v^fd*! 


( 10 ) 


Also, by expressing this quantity in terms of a 3+1 decomposition of a symmetric second rank tensor leads 
to the same expression as in Ref. |14] . that is 

( 11 ) 

where £ is the internal energy density, p the hydrostatic pressure, q°‘ the heat flux and 11"^ Navier’s viscous 
tensor. Also, is the so-called spatial projector which is given by 

The force term on the right hand side of Eq. ([9]) can be expressed as the four-divergence of the electromagnetic 
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stress tensor PCS] such that the energy-momentum conservation equation for the single charged fluid 
can be written as 


A. 

dx^ 


(T^“ + = 0 


(13) 


Equation (IT^ contains both the momentum balance equation and the total energy balance. In order to write 
an expression for the internal energy evolution one projects Eq. in the direction of the hydrodynamic 
velocity namely 


d 


+ M«^) = 0 


which yields 


dW ,1 dq'' 


(14) 


(15) 


dx’^ ^ dx'^ 

where the dot indicates a fluid’s proper time derivative, i.e. Notice that the electromagnetic 


field term vanishes since 


= 0 


(16) 


where we have used that is antisymmetric. For the momentum balance one obtains, from Eq. (lT3l) . 




(17) 


Equations m, (USD and fll7p constitute the set of transport equations for a relativistic single component 
inviscid fluid in the presence of an electromagnetic field . A constitutive equation that expresses the heat 
flux in terms of the forces is required in order to close the set such that it can appropriately describe the 
corresponding dynamics. The Chapman-Enskog procedure will be used in the following sections in order to 
obtain such a relation. 


3 The electromagnetic contribution to 

In this section, the Chapman-Enskog method is used in order to establish a hrst order in the gradients 
correction to the equilibrium distribution function using Eq. CDi. For the system here considered, all gradients 
of the hydrodynamic velocity are negelcted since the heat flux, being a vector flux, will not be coupled to 
second rank tensor forces. 

For the system in hand, a simple charged gas in the presence of an electromagnetic field, the acceleration 
acting on the molecules that appears in the second term of Eq. (??) is calculated using the electromagnetic 







field tensor which, in cartesian coordinates reads 
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Also, for the sake of simplicity, a relaxation time model will be introdnced on the right hand side of Boltzmann’s 
eqnation by assnming that all details of collisions can be inclnded in one parameter r snch that 


jun = 


f - /<“> 


r 


(19) 


This simple model has been extensively nsed in order to assess the strnctnre of the flnxes-forces relations for 
single component flnids. For a more precise valne of the transport coefficients present in snch relations, the nse 
of the complete kernel is reqnired. Introdncing Eqs. (IT^ and (fT8|) in Eq. ([1]) together with Chapman-Enskog’s 
hypothesis, namely 

/ = ( 20 ) 


one obtains for the non-eqnilibrinm correction to the distribntion fnnction 




( 21 ) 


where 


/'”> = 


n 


-e 


4:7r(^zK2 (^) 

is the local eqnilibrinm solntion j2] which leads to the expression [8] 


( 22 ) 


/W = -T {,.“/?> 

L ’ \mc^z 




(23) 


In Eq. fl2^ Q) is the n — th order modified Bessel fnnction of the second kind. On the other hand, 
for the first term on the right hand side of Eq. (|2T1) . the local eqnilibrinm Enler eqnations need to be 
introdnced in order to write the time derivatives of the state variables, that arise from the fnnctional hypothesis 
/ {x'^, v'") = / (n*^, n [x'^) , T (x'^), U {x'^)), in terms of the forces. In order to separate time and space derivatives 









in such a term we introduce the following decomposition for an arbitrary reference frame 


= - 




P \ ^2 
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( 24 ) 


such that 
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For the hrst term in Eq. fl3^ one has 


(25) 
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where Q (^) = K 3 (i) !K 2 (^) and for the second one 


m“/,2> = /'“> 
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zc^ 
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K"Ta + 

zc^ 


(26) 


(27) 


As mentioned above, the proper time derivatives are eliminated intriducing Euler’s equations which are 
obtained by neglecting dissipative terms in the transport equations (IH]), (1151) and flTTl) . which are given by 


,,,, dn dhi^ 

W -— = -n- 


dx^ 


dxf^ 


(28) 




= 


P u'' 


(29) 

(30) 


where the specihc heat at constant number density Cn = (§f)n been introduced. Introducing Eqs. 
in Eq. (l39l) yields the following expression 


y«/(o) = /(o) 
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(31) 


Notice that the force term in Eq. (|3T1) will compete with the second term on the right hand side of Eq. 
Isolating these two electromagnetic contributions to the deviation from equilibrium of the distribution 
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function due to them, which we call f^lf, 


we obtain 


f(i) _ 
JEM ~ 



(32) 


where we have used p = nmQ (i). Equation (132]) is the main result of this section. It consitutes the 
generalization of the result obtained in Ref. j8] where the fact that a purely electrostatic field can drive a 
single charged gas out of equilibrium in a relativistic scenario. In the present work, the generalization is 
twofold in the sense that it includes a magnetic held and is consequently expressed in an arbitrary reference 
frame. This expression will be used in the next section in order to establish the heat hux due to this effect. 


4 Calculation of the heat flux 


As extensively discussed in Refs. na and jT9], heat is dehned as the transport of internal energy arising from 
the chaotic motion of the molecules which is in principle a quantity characteristic of the comoving frame in 
which all mechanical effects are absent. Also, as shown in Ref. mi, it can be constructied as a tensor quantity 
and can thus be calculated in any arbitrary frame. The expression 


Q" = Ajr’ / k-’phmd'K 


(kp 


(33) 


is here used in order to assess the contribution to thermal dissipation due to the electromagnetic held mi 
in an arbitrary reference frame in which the hydrodynamic four velocity of the huid is given by In Eq. 


is a Lorentz boost with velocity u and K'^ = 7 (fc) 


k, c 


is the chaotic velocity, that is the particles’ 


velocities measured by an observer comoving with the huid’s volume element 


In order to calculate the contribution to the integral in Eq. (|33l) arising from h needs to be written 
in terms of the chaotic velocity. In order to accomplish such task, we use the invariant expression 


v^Ub 


-l(k) 


(34) 


and also write the molecular velocity as 


= 'y(k)jO,i,uk‘' 


(35) 
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Introducing these expressions in Eq. fl32|) one obtains 


f(i) _ 
Jem — 




zmc 


6(M 


-1 


(36) 


Thus, the electromagnetic contribution to the heat flux is given by 


Qe« = j - l| 7,%d*A- 

which can also be written as 

QU = (^) 

Z 

where the integral X"'*' is given by 


(37) 


(38) 


r’ (z) = /<")rV - l|7,t,\/7?,,-l<i7 (39) 

Equation (l38|) contains the main result of this work: the heat flux for a single component charged fluid in 
the special relativistic regime features a contribution due to the electromagnetic force This effect 

is purely relativistic and generalices the result obtained in Ref. [ 8 ] where the particular case of a purely 
electrostatic held in a comoving frame was addressed. It is important to emphasize that the dissipative flux 
in Eq. (j38l) is a tensor quantity. 

In order to assess the magnitude of the corresponding transport coefficient, the values for the integral 
jafe need to be obtained. This can be readily done yielding 


[z) = 0, for a 7 ^ 6 , a, 6 = 1, 2, 3,4 
X““ (z) = nz(? ( 5^ + ^ j for a = 1, 2, 3 




Also notice that the 7 = a = 4 term in Eq. (1381) vanishes since £^4 cx and thus one can write 


Q%m- ''™(iT)5 




(40) 


with a,b = 1,2,3. The non-relativistic thermal conductivity in the relaxation time approximation k^r = 


Q 











has been here introduced. 

2 m 

The components of the heat flux in Eq. flTOj) can be explicitly obtained in an arbitrary frame where 
W’' = [u, c], which yieds: 


Qem - -'t-VH 5 I 5 + - — 




E + u X B, E ■ 


u 


from which one can extract three relevant limits. Firstly notice that in the comoving frame where U'' = 
and 7 u = 1 one obtains 


(41) 


0 , c 




E,0 


which is precisely the result obtained in Ref. [8] for the case of a purely electrostatic field in the fluid’s 
comoving frame. The other two limits pertain the non-relativistic case. For a fluid with a high value of but 
such that the hydrodynamic velocity is considerably lower than the speed of light, ~ 1 and 
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(42) 


which is physically meaningfull since the thermodynamic force in this limit is the electromagnetic force for 
the spatial components and corresponds to the mechanical electrostatic dissipation on the fourth component. 
Finally, the striclty non-relativistic limit in which also the temperature is low such that z -C 1 one has 
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EM 


rxj 


— I^NR 



E + u X B 



finally showing that this effect is strictly relativistic. 


(43) 


5 Final remarks 

This paper has been devoted to the extension of the Benedicks-type effect recently identified in the context 
of special relativity to the magnetic field case. Although the calculation may seem to be a straightforward 
exercise considering the inclusion of the magnetic field components in Faraday’s field tensor when calculating 
there are several features that must be highlighted. The fact that B appears explicitly in implies 
a non-zero contribution of the magnetic field to the entropy production present in a single component fluid. 
In contrast, in the non-relativistic case the magnetic field does not contribute to the entropy production at 
the Navier-Stokes level |18| . More over the hydrodynamic equations now include new terms in which the field 



















is now identified as a thermodynamic force. This is rather new, since the only dissipative effects related to 
magnetic fields previously identified at first order in the Knudsen parameter (for simple fluids) were those 
associated to the values of the viscosity coefficients [IH1[I21E3]. 

Magnetic fields are well known to be ubiquous in the universe. Although the corresponding astrophysical 
mean free times are quite large compared to the ordinary scales, the characteristic times for most processes 
are such that dissipative effects must be taken into account [Spitzer]. It is in this type of scenarios in which 
irreversible thermodynamics becomes irrelevant. Future work will be devoted the study of hydromagnetic 
instabilities including the effects presented in this work. 
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